Role of contact formation process in transport properties of molecular junctions: 
conductance of Au/BDT/Au molecular wires 
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We report theoretical investigations on the role of contact formation process and its resulting 
structures to quantum transport in molecular wires and show that these processes critically control 
charge conduction. It was found, for Au(lll)/l,4-benzenedithiol(BDT)/Au(lll) junctions, the 
hydrogen atom in the thiol groups is energetically non-dissociative after the contact formation. The 
calculated conductances and junction breakdown forces of H- non-dissociative Au/BDT/Au devices 
are consistent with the experimental values, while the H-dissociated devices give conductances more 
than an order of magnitude larger. The results can be well understood by examining the scattering 
states that traverse the contacts. This work clearly demonstrates that the contact formation process 
must be carefully understood in order to correctly capture quantum transport properties of molecular 
nanoelect ronics . 

PACS numbers: 73.63.-b, 73.63.Rt 85.35.-p, 73.40.Ns 



In a very recent paper, Song et al.[l\ reported experi- 
mental fabrications of three-terminal single molecule field 
effect transistors. The successful fabrication and charac- 
terization of such a molecular device can be considered an 
important milestone of nanoelectronics. In Song et a/.[l[ 
experiment, a gold wire was broken by electro- migration 
to produce a nano-meter gap in the wire. A molecule 
such as l,4-benzenedithiol(BDT) may bridge the gap and 
form an Au/BDT/Au single molecule transport junc- 
tion. The equilibrium conductances of 0.01 ~ O.OlSGo 
(Go = 2e^/h) were reportedflj giving an average value 
of (0.0132 ± 0.0021)Go. The Au/BDT/Au device was 
actually subjected to extensive studies in the past, for 
instance Xiao et at. reported a mean conductance of 
0.011 Go by statistically measuring several thousand 
Au/BDT/Au junctions formed by BDT bridging the gap 
between a scanning tunneling microscope (STM) and a 
Au surface [H. Their result has also been reproduced by 
other experiments [3i] . The consistency of measured trans- 
port properties on devices fabricated by totally different 
methods suggests a degree of structural-function control 
at the single molecule level. 

Despite the achievements, one realizes that the atomic 
structures of the fabricated devices were still unknown. 
In particular, the most important structural information 
- the metal- molecule contacts, is at best ambiguous for 
essentially all single molecule transport junctions inves- 
tigated in literature. In general, a most important sci- 
ence issue concerning nano-systems is the relationship be- 
tween structure and function. The experimental conver- 
gence of transport data for the Au/BDT/Au device pro- 
vides a timely opportunity to shed considerable light on 
the structure- function issue of molecular nanoelectronics. 
It is the purpose of this Letter to report our first princi- 
ples theoretical investigation on how the contact forma- 
tion process can critically affect quantum transport prop- 



erties of molecular junctions. In all theoretical analysis of 
molecular devices, one assumes an initial contact struc- 
ture between the molecule and the metal electrodes — 
guided by intuition or by experiments, and then relaxes 
the structure. However, experimentally when a molecule 
is brought to contact the metal leads, a contact forma- 
tion process occurs where chemical reactions may give 
rise to dissociation or formation of atomic groups from 
the original molecule. Such a process is likely lost when 
an initial atomic configuration is assumed without care- 
fully considering it from ah initio point of view. As a 
result, the formation process has not been subjected to 
systematic investigations so far and, as we show below, it 
is a crucial effect that controls the interface transparency 
to charge flow. 

We use the Au(lll)/l,4-benzenedithiol(BDT) / 
Au(lll) as a prototypical system for our investigation. 
As discussed above, the experimental conductance value 
is O.Oll-O.O15Go0,l3 for this system. On the theoretical 
sidejJ-Q? conductance obtained from quantitative analy- 
sis have not been able to reach a consensus. In particular, 
density functional theory (DFT) based first principles 
methods within local density approximation (LDA) or 
generalized gradient approximation (GGA) have mostly 
produced conductance values considerably larger — by 
more than one or even two orders of magnitude, than the 
experimental value of Xiao et a/.Q and Song et al\l\. In 
this Letter, we shall show that such a qualitative dis- 
crepancy can be well understood by investigating the 
role of the device contact formation process. By exten- 
sive first principles calculations, it was found that a BDT 
molecule prefers to attach to ad-atoms when it is bonded 
to Au(lll) and, much more importantly, the hydrogen 
atoms on the thiol groups of a BDT do not dissociate 
away after the formation of the device contacts. The cal- 
culated conductance of H- non-dissociated Au/BDT/Au 



junctions as well as the junction breakdown force are 
consistent with experimentally measured values. On the 
other hand, all the H-dissociative contact configurations 
produce conductances at least one order of magnitude 
greater than the experimental value of Refs.[l|, |2|, while 
the junction breakdown force is also significantly larger 
than the experimental results. The quantum transport 
results can be well understood by analyzing the scatter- 
ing states that traversing the metal- molecule contacts. 
Bonding structure of Au/BDT interface. De- 
spite the importance of bonding structure of the metal- 
molecule interface to transport [s', detailed bonding 
geometry of Au/BDT/ Au junctions have not been well 
established. In different systems such as the popular 
Au/alkanethiol/Au junction, the conductance is dom- 
inated by electronic tunneling through the localized 
cr— bonds of the alkaneflo), namely the huge resistance 
(typically tens of mega-Ohms) of Au/alkanethiol/Au is 
dominated by the length of alkane molecule and not 
by the Au/alkanthiol interface. In Au/BDT/ Au, the 
benzene ring in BDT consists of a non-local tt— bond 
system which is rather transparent to electron conduc- 
tion hence the conductance should be much more sen- 
sitive to the bonding structures at the Au/BDT inter- 
face. A very careful calculation on the interface is thus 
necessary. To this end, we use a standard DFT-PAW 
method with the Perdew-Burke-Ernzerhof 96 functional 
(GGA-PBE) as implemented in the electronic structure 
package VASP[11] to determine the atomic structure [12]. 
As suggested by a recent STM experiment [HI , thiol 
molecules prefer to attach Au(lll) surface through Au 
ad-atoms. We have therefore calculated BDT absorption 
on Au(lll) with and without Au ad- atoms. It was found 
that if there are ad-atoms, the structures with BDT at- 
taching to the surface via ad-atoms are more stable than 
their direct adsorption to the surface: the difference of 
binding energies is at least 0.4 eV per molecule. As a 
result we shall focus on a series of representative atomic 
configurations with BDT molecules absorbed on Au(lll) 
surface via Au ad-atoms. Figure [Ha,b) plots two ini- 
tial structures of BDT attached to Au(lll) via an ad- 
atom prepared for the structural relaxations, where the 
molecule is parallel or perpendicular to the Au(lll). Ac- 
cording to the C3 rotational symmetry of Au(lll), an- 
other two sets of these structures with a rotation angle of 
—30^ or 30^ were also considered. Together these should 
likely cover most of the initial configurations. For the 
contact formation process, we consider three cases for 
the dissociative configuration: the dissociated H atom 
attaches to the ad-atom, to the surface, or escapes to vac- 
uum to form an H2 molecule, as shown in Fig. [I](c). The 
calculated total energies are summarized in Tab. I. It is 
striking to find that configurations with non-dissociative 
S-H bonds, (first row in Tab. I) are always energetically 
more stable than all dissociative structures by at least 
0.2 eV per BDT throughout all systems investigated. In- 




FIG. 1: (Color online) Top views of examples for (a) paral- 
lel and (b) perpendicular configurations, (c) Four positions 
of H have been considered for each configuration at several 
orientations: (A) is for non-dissociated H; (B,C,D) are for 
dissociated H which attaches to the (B) ad-atom , (C) the 
surface , or (D) escapes into vacuum forming an H2 molecule. 



parallel perpendicular 

0° 30° 60° 0° 30° 60° 

H-non-dissociative 0.00 0.09 0.08 0.06 0.06 0.06 

H-dissociated to ad-atom 0.28 0.25 0.61 0.53 0.54 0.42 

H-dissociated to surface 0.41 0.40 0.42 0.40 0.41 0.41 

H-dissociated to vaccum 0.24 0.25 0.26 0.25 0.25 0.24 

TABLE I: The difference of total energies (units eV) compared 
to the the most stable structure (the parallel 0°) for typical 
configurations of Au/BDT interface. For all situations, the 
H-non-dissociative structures have lower energies. 



deed, the non-dissociative S-H bond in thiol group ab- 
sorbed on a perfect Au(lll) surface was experimentally 
observed by Yates et al. |l4| and more recently confirmed 
by ab initio calculations [15]. Furthermore, our calcu- 
lated breakdown force of the junction (by elongation) is 
from 1.1 nN to 1.6 nN[l6|, fairly consistent to experimen- 
tal measurements of molecular junctions with the same 
thiol linker [17]. These results strongly suggest that an 
H-non-dissociative structure should provide a more real- 
istic model in terms of transport modeling of the BDT 
devices. 

Bonding structure in Au/BDT/ Au device. The 

parallel configurations were adopted to build the two- 
probe model of the junction, since they have lower ener- 
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FIG. 2: (Color online) Atomic models of the non-dissociated 
BDT linked by ad-atoms (a) without and (b) with stretching; 
and the H dissociated BDT connected to the surface via (c) 
ad-atoms or (d) Au atomic chains of various lengths. 
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FIG. 3: (Color online) Transmission of H-dissociative (solid 
blue) and H-non-dissociative (dashed red) models versus en- 
ergy without junction stretching. The cross denotes the ex- 
perimental value 0.0 11 Go |2] . Inset: conductance (Co) versus 
junction length (A) under mechanical stretching (Fig. ^jp) for 
H-non-dissociative model. 



gies than the perpendicular ones, as shown in Fig. [2fa). 
In our transport calculations (below), a two-probe model 
of Au/BDT/Au device is comprised of three parts [Tsj: 
the scattering region plus the left /right leads. The scat- 
tering region includes a BDT molecule bonded to four 
layers of Au(lll) atoms via an ad- atom on either side 
(atoms inside the black rectangular box in Figl2j). The 
leads are bulk-like Au(lll) layers extending on either 
side of the scattering region to z = ±00 where z is 
the direction of current flow. Periodic boundary con- 
dition is applied in x — y directions. In order to simulate 
the pulling process of the junctions in the experiment of 
Xiao et a/. [2], we have constructed a series of junctions 
with different distances L between the two surfaces of 
the left/right electrodes (see Ref.[16]). For comparison 
purposes, we have also investigated the widely assumed 
H-dissociated model. When the S-H bond is dissociated^ 
the Au-S bond shows a stronger strength than the Au- 
Au bond in our calculation. Therefore, in the experi- 
ment of Xiao et al. where a STM tip was repeatedly re- 
tracted from a Au surface, it is quite possible to drag an 
Au atomic wire out of the surface through the attached 
BDT molecule [2]. To cover this possibility, for the H- 
dissociative structure, we have calculated a few devices 
where the H-dissociated BDT is linked to the surface 
via gold atomic chains with various lengths (Fig. [2li). 
The geometries of the scattering region of all instances of 
Au/BDT/Au are again fully optimized [12] using VASP. 
It was found that the H- non-dissociated junctions al- 
ways have lower total energies than the corresponding 
H-dissociated cases: because there are two Au/BDT in- 
terfaces now, the energy difference roughly doubles those 
listed in Tab. 1. 

Transport properties. Our transport analysis is 
based on carrying out DFT analysis within the Keldysh 
non-equilibrium Green's function (NEGF) formalism [18|. 
The basic idea of NEGF-DFT is to self-consistently cal- 
culate the Hamiltonian of the device by DFT and deter- 




FIG. 4: (Color online) Scattering states of: (a) H-dissociative; 
(b) H-non-dissociative models. Insets are the HOMO and 
LUMO of the corresponding molecule 1,4-benzenedithiolate 
and 1,4-benzenedithiol. 



mine the non-equilibrium quantum statistical properties 
of the device operation by NEGF. For more details we 
refer interested readers to the original literature [18]. In 
the NEGF-DFT self-consistent calculation of the density 
matrix and Hamiltonian, we use double-^ plus polariza- 
tion(DZP) linear combination of atomic orbital (LCAO) 
basis sets for all the atoms, GGA-PBE for the exchange- 
correlation potential, and define atomic core potentials 
using standard norm conserving pseudopotentials. 

Figure [3] shows the calculated transmission(T) spectra 
of Au/BDT/Au devices versus energy in an energy 
range of —3.0 eV to 3.0 eV. For the H-non-dissociated 
model, the equilibrium conductance (the value of T at 
Fermi level) was found to be 0.065Go for a junction at 
its equilibrium junction length (L = 13.8 A, see Ref.[16]). 
This conductance decreases when stretching the junc- 
tion, and reaches 0.022Go for L ^ 15.7 A, at which 
the junction starts its mechanical breakdown (see EPAPS 
file[16] for details). The inset of Fig. [3]plots conductance 
changes versus L. This range of conductance is within 
the same order of magnitude as the experimental value 
of 0.011Go[2] to 0.0132GopL]. In fact, it is most likely 
that the stretched junction is closer to the experimen- 
tal reality since during the measurements [2], a STM tip 
was retracting away from the Au surface thereby stretch- 
ing the molecular wires. Hence the calculated conduc- 
tance of O.O22G0 just before junction breakdown is quite 
reasonable to theoretically describe the transport. In 
comparison, the conductance of all H-dissociated models 
(Fig. [2][c) and (d)) with or without stretching are much 
higher, in a range of 0.38 — 0.86Go- These high values are 
consistent with previous ah initio calculations [sl-Q which 
assumed H-dissociated models. Therefore, the H-non- 
dissociated Au/BDT/Au junctions which are energeti- 
cally more stable, have conductance values much closer 
to the experimental data — by one order of magnitude, 
than the H-dissociated models. 

The transport results clearly suggest that the hy- 
bridization of electronic states from Au electrodes 
and the molecule is significantly different for H-non- 
dissociative and dissociative models. To understand why 
this happens, we have analyzed the scattering states 
around the Fermi level for both models, as shown in Fig. 
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m When the H is dissociated from a S-H group, an elec- 
tron of the S atom becomes unpaired which has an over- 
whelming tendency to attract an additional electron to 
make a pair. The additional electron is most likely con- 
tributed by s-electrons of Au leads, resulting in a trans- 
fer of charge from leads to the molecule. The transferred 
electron dopes into the lowest unoccupied molecular or- 
bital (LUMO) of the molecule and pushes down the 5- 
LUMO bonding statejust below the Fermi level, as found 
in a similar system[l9]. The 5-LUMO bonding states are 
expected to be very delocalized, since it was composed 
of a delocalized LUMO and a delocalized s state. We 
thus plot the scattering states around the Fermi level, as 
shown in Fig. IH^a). It was found that the conductance 
around the Fermi level is indeed dominated by a delocal- 
ized 5-LUMO state as expected. The plot shows the scat- 
tering states nicely passing through the junction, giving 
rise to a high conductance value. On the other hand, for 
the H-non-dissociative model, the above charge transfer 
can hardly happen since all electrons are paired already, 
hence the hybridized 5-LUMO state around the Fermi 
level disappears. The bonding picture then switches to 
that of a lone-pair of the S atom donating to the partially 
unfilled 5-band of Au leads. The lone-pair is rather lo- 
calized, like a a-type orbital. The probable hybridization 
between the lone-pair and Au leads is therefore some- 
what localized, which results in a tunneling mechanism 
for electrons going through the junction at low bias: a 
much smaller conductance is therefore expected. Indeed, 
as shown in the plot of scattering states in Fig. IH^b), very 
few incoming scattering states can pass through the junc- 
tion. In this case, the conductance is mainly contributed 
by the HOMO of 1,4-benzenedithiol. 

In summary, on the structural side our investigations 
reveal, surprisingly, that the H atoms in the thiol group 
of the BDT are non-dissociative after the Au/BDT/Au 
transport junctions are formed. On the functional side, 
the introduction of non-dissociative H atoms blocks 
charge transfer doping to the BDT from the Au elec- 
trodes, effectively induces an extra potential barrier that 
considerably reduces the electron transparency of the 
Au/BDT interface. Since thiol molecules provide per- 
haps the most popular binding linkers in experiments, 
these findings shed considerable light on charge conduc- 
tion properties at the single molecule level. For the H- 
non-dissociative structures, we predict a conductance of 
O.O22G0 for stretched Au/BDT/Au junctions before me- 
chanic breakdown, ^hich is in good consistency to the 
measured value [l, 2]. In comparison, all H-dissociated 
junctions produce conductance values more than one or- 
der of magnitude greater. Furthermore, the junction 
breakdown force for the H-non-dissociative structure is 
also consistent with the measured value [Tt), and that 
for H-dissociative system is significantly higher. Our in- 
vestigation clearly reveals how structure formation could 
drastically infiuence transport properties. Finally, we be- 



lieve that our findings resolved a long-standing theory- 
experiment discrepancy on the conductance of the most 
widely studied molecular transport junction. 
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